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We present a polynomial time algorithm to compute any fixed number of the highest 
coefficients of the Ehrhart quasi-polynomial of a rational simplex. Previously such 
algorithms were known for integer simplices and for rational polytopes of a fixed 
dimension. The algorithm is based on the formula relating the fcth coefficient of 
the Ehrhart quasi-polynomial of a rational polytope to volumes of sections of the 
polytope by affine lattice subspaces parallel to fc-dimensional faces of the polytope. 
We discuss possible extensions and open questions. 



1. Introduction and main results 

Let P C M. d be a rational polytope, that is, a polytope with rational vertices. 
Let t e N be a positive integer such that the dilated polytope 

tP = {tx: x E P\ 



has integer vertices. As is known, see, for example, Section 4.6 of [St97], there exist 
functions ej(P; •) : N — ► Q, i = 0, . . . , d, such that 

ei(P;n + t) = e»(P;n) for all neN 

and 

d 

iP nZ d \ = y^ei(P;n)n* for all neN. 

i=0 



ru 



The function on the right hand side is called the Ehrhart quasi-polynomial of P. It 
is clear that if dimP = d then ed(P; n) = vol P. In this paper, we are interested in 
the computational complexity of the coefficients ei(P;n). 
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If the dimension d is fixed in advance, the values of e$(P;n) for any given P, n, 
and i can be computed in polynomial time by interpolation, as implied by a poly- 
nomial time algorithm to count integer points in a polyhedron of a fixed dimension 
[B94a], [BP99]. 

If the dimension d is allowed to vary, it is an NP-hard problem to check whether 
P H Z d 7^ 0, let alone to count integer points in P. This is true even when P is a 
rational simplex, as exemplified by the knapsack problem see, for example, Section 
16.6 of [Sc86]. If the polytope P is integral then the coefficients ei(P;n) = ei{P) 
do not depend on n. In that case, for any k fixed in advance, computation of 
the Ehrhart coefficient ed-k(P) reduces in polynomial time to computation of the 
volumes of the (d — /c)-dimensional faces of P [B94b]. The algorithm is based 
on efficient formulas relating e^_fc(P), volumes of the (d — /c)-dimensional faces, 
and cones of feasible directions at those faces, see [Mo93], [BP99], and [PT04]. In 
particular, if P = A is an integer simplex, there is a polynomial time algorithm for 
computing ed_fc(A) as long as k fixed in advance. 

In this paper, we extend the last result to rational simplices. 

• Let us fix an integer k > 0. The paper presents a polynomial time algorithm, 
which, given an integer d > k, a rational simplex A C M d , and a positive integer n, 
computes the value of ed-fc(A;n). 

We present the algorithm in Section 7 and discuss its possible extensions in 
Section 8. 

In contrast to the case of an integral polytope, for a general rational polytope P 
computation of ej(P; n) cannot be reduced to computation of the volumes of faces 
and some functionals of the "angles" (cones of feasible direction) at the faces. A 
general result of McMullen [Mc78], see also [MS83] and [Mc93], asserts that the 
contribution of the i-dimensional face F of a rational polytope P to the coefficient 
ei{P;n) is a function of the volume of F, the cone of feasible directions of P at F, 
and the translation class of the affine hull aff (F) of F modulo Z d . 

Our algorithm is based on a new structural result, Theorem 1.3 below, relating 
the coefficient ed-k(P',n) to volumes of sections of P by affine lattice subspaces 
parallel to faces F of P with dimF > d — k. Theorem 1.3 may be of interest in its 
own right. 

(1.1) Valuations and polytopes. Let V be a (i-dimensional real vector space 
and let A C V be a lattice, that is, a discrete additive subgroup which spans V. 
A polytope P C V is called a A-polytope or a lattice polytope if the vertices of P 
belong to A. A polytope P C V is called A-rational or just rational if tP is a lattice 
polytope for some positive integer t. 

For a set A C V, let [A] : V — > R be the indicator of A: 
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A complex-valued function v on rational polytopes P C V is called a valuation if 
it preserves linear relations among indicators of rational polytopes: 

iei iei 

where Pi C V is a finite family of rational polytopes and cti are rational numbers. 
We consider only A-valuations or lattice valuations v that satisfy 

v{P + u) = u(P) for all u G A, 

see [MS83] and [Mc93]. 

A general result of McMullen [Mc78] states that if v is a lattice valuation, P G V 
is a rational polytope, and t G N is a number such that tP is a lattice polytope 
then there exist functions Vi{P\ •) : N — >■ C, i = 0, . . . , d, such that 

d 

u(nP) = ^2 v ii p 'i n ) ni for a11 n eN 

i=0 

and 

Vi(P;n + t) = Ui(P;n) for all n G N. 

Clearly, if we compute v{mP) for m = n, n + t, n + 2t, . . . , n + dt, we can obtain 
Vi(P;n) by interpolation. 

We are interested in the counting valuation E, where V = M d , A = Z d , and 

E{P) = \PHZ d \ 

is the number of lattice points in P. 

The idea of the algorithm is to replace valuation E by some other valuation, so 
that the coefficients ed(P',n), . . . ,ed_A;(P;n) remain intact, but the new valuation 
can be computed in polynomial time on any given rational simplex A, so that the 
desired coefficient ed-fc(A;n) can be obtained by interpolation. 

(1.2) Valuations El. Let L C M. d be a lattice subspace, that is, a subspace 
spanned by the points LflZ d . Suppose that dimL = k and let pr : IR d — > L be the 
orthogonal projection onto L. Let P C M. d be a rational polytope, let Q = pr(P), 
Q C L, be its projection, and let A = pr (Z d ) . Since L is a lattice subspace, Ac I 
is a lattice. 

Let L 1 - be the orthogonal complement of L. Then C l d is a lattice sub- 
space. We introduce the volume form vold-k on L which differs from the volume 
form inherited from R d by a scaling factor chosen so that the determinant of the 
lattice Z d n L is 1. Consequently, the same volume form vo\d-k is carried by all 
translations x + lA, x G M d . 
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We consider the following quantity 



E L (P)=J2vol d - k {Pn{m + L ± ))= vol d _ fc (PR (m + L^)) 

m£A meQnA 

(clearly, for m Q the corresponding terms are 0). 

In words: we take all lattice translates of L- 1 , select those that intersect P and 
add the volumes of the intersections. 

Clearly, E L is a lattice valuation, so 



E L (nP) =^e,(P,L;n)n l 



i=0 

for some periodic functions ej(P, L; •). If tP is an integer polytope for some teN 
then 

ei(P,L;n + t) = ei(P,L;n) for all n£l 

and i = 0, . . . , d. 

Note that if L = {0} then E L {P) = volP and if L = R d then E L (P) = |PnZ d |, 
so the valuations El interpolate between the volume and the number of lattice 
points as dimL grows. 

We prove that ed-k(P',n) can be represented as a linear combination of 
ed-h(P, L; n) for some lattice subspaces L with dimL < k. 

(1.3) Theorem. Let us fix an integer k > 0. Let P C M d be a full- dimensional ra- 
tional polytope and let t be a positive integer such that tP is an integer polytope. For 
a (d — k) -dimensional face F of P let \m(F) C M. d be the (d — k)- dimensional sub- 
space parallel to the affine hull aff(F) of F and let L F = (linF)" 1 be its orthogonal 
complement, so L F C M d is a k-dimensional lattice subspace. 

Let C be a finite collection of lattice subspaces which contains the subspaces L F 
for all (d — k)- dimensional faces F of P and is closed under intersections. For 
L G C let u(L) be integer numbers such that the identity 



Lee . 



Lec 



holds for the indicator functions of the subspaces from C 
Let us define 

v(nP) = v(L)E L (nP) for 
Lec 

Then there exist functions Vi(P; •) : N — > Q, i = 0, . . . , d, such that 

d 

(1.3.1) u(nP) = ^2ui(P; ri)n % for all tiGff, 

i=0 
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(1.3.2) 



Vi(P\n + t) = h>i(P;n) for all tiGN, 



and 

(1.3.3) ed-i(P;n) = Ud-i(P;n) for all tiGN and i = 0,...,k. 

We prove Theorem 1.3 in Section 4 after some preparations in Sections 2 and 3. 

The advantage of working with valuations El is that they are more amenable 
to computations. 

• Let us fix an integer k > 0. We present a polynomial time algorithm, which, 
given an integer d > k, a (i-dimensional rational simplex A C M d , and a lattice 
subspace L C M d such that dimL < k, computes El(A). 

We present the algorithm in Section 6 after some preparations in Section 5. 

(1.4) The main ingredient of the algorithm to compute ed_fc(A;n). 

Theorem 1.3 allows us to reduce computation of ed-k{A;n) to that of -Ex(A), 
where L C M. d is a lattice subspace and dimL < k. Let us choose a particular 
lattice subspace L with dimL = j < k. 

If P = A is a simplex, then the description of the orthogonal projection Q = 
pr(A) onto L can be computed in polynomial time. Moreover, one can compute in 
polynomial time a decomposition of Q into a union of non-intersecting polyhedral 
pieces Q i: such that voLj-j (pr~ 1 (x)) is a polynomial on each piece Qi. Thus 
computing of E L (A) reduces to computing of the sum 

meQiDA 

where <f> is a polynomial with deg^> = d — j, Qi C L is a poly tope with dimQi = 
j < k and A C L is a lattice. The sum is computed by applying the technique of 
"short rational functions" for lattice points in polytopes of a fixed dimension, cf. 
[BW03], [BP99], and [D+04]. 

The algorithm for computing the sum of a polynomial over integer points in a 
polytope is discussed in Section 5. 

2. The Fourier expansions of E and E l 

Let V be a (i-dimensional real vector space with the scalar product (•, •) and the 
corresponding Euclidean norm || ■ ||. Let A C V be a lattice and let A* C V be the 
dual or the reciprocal lattice 



A* = |x G V : (x, y) e Z for all y E A 
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For r > 0, we introduce the theta function 



9a{x, t) =T d l 2 exp {— 7rr||x — m|| 2 } 

m£A 

= (detA) _1 ^ exp{-7r||/|| 2 /r + 2ttz(/,x)} , where x G V. 
leA* 

The last inequality is the reciprocity relation for theta series (essentially, the Poisson 
summation formula), see, for example, Section 69 of [Be61]. 

For a polytope P, let int P denote the relative interior of P and let dP = P\int P 
be the boundary of P. 

(2.1) Lemma. Let P C V be a full- dimensional polytope such that dP D A = 0. 
Then 



PnA 




Proof. As is known (cf., for example, Section B.5 of [La02]), as r — > +oo, the 
function 0\{x,t) converges in the sense of distributions to the sum of the delta- 
functions concentrated at the points m G A. Therefore, for every smooth function 
(j) : M. d — > R with a compact support, we have 

(2.1.1) lim / </>(x)6a(x,t) dx = 0(m). 

Since <9Pfl A = 0, we can replace <p by the indicator function [P] in (2.1.1). □ 

Remark. If dP n A ^ 0, the limit still exists but then it counts every lattice point 
m G dP with the weight equal to the "solid angle" of m at P, since every term 
exp {— txt\\x — m|| 2 } is spherically symmetric about m. This connection between 
the solid angle valuation and the theta function was described by the author in the 
unpublished paper [Ba92] (the paper is very different from paper [B94b] which has 
the same title) and independently discovered by Diaz and Robins [DR94] . Diaz and 
Robins used a similar approach based on Fourier analysis to express coefficients of 
the Ehrhart polynomial of an integer polytope in terms of cotangent sums [DR97] . 
Banaszczyk [B93a] obtained asymptotically optimal bounds in transference theo- 
rems for lattices by using a similar approach with theta functions, with the polytope 
P replaced by a Euclidean ball. 

The formula of Lemma 2.1 can be considered as the Fourier expansion of the 
counting valuation. 

We need a similar result for valuation El defined in Section 1.2. 
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(2.2) Lemma. Let P C M. d be a full- dimensional polytope and let L C R d be a 
lattice subspace with dimL = k. Let pr : M d — > L be the orthogonal projection 
onto L, let Q = pr(P), and let A = pr(Z d ), so A C L is a lattice in L. Suppose 
that dQ n A = 0. 
Then 

E L (P)= lim V exp{-7r||/|| 2 /r} / exp{2yrz(/, x)\ dx. 
ieLnz d ^ 

Proof. We observe that LflZ d = A*. For a vector x G M. d , let xl be the orthogonal 
projection of x onto L. Applying the reciprocity relation for theta functions in L, 
we write 

exp{-7r||/|| 2 /r + 27rz(/,x)} 

ieLnz d 

= exp{-7r||/|| 2 /r + 27rz(/,x L )} 

ieLnz d 

= (det A)r fc/2 ^ exp {--kt\\x l - m\\ 2 } . 

As is known (cf., for example, Section B.5 of [La02]), as r — ► +00, the function 

g T (x) = T k l 2 exp {— tvt\\xl — m|| 2 } 

m£A 

converges in the sense of distributions to the sum of the delta- functions concentrated 
on the subspaces m + L 1 - (this is the set of points where xl = m) for m E A. 

Therefore, for every smooth function <fi : R d — >■ H. with a compact support, we 
have 

(2.2.1) lim / 4>{x)g T {x) dx = V ] / <f>(x) d L ±x, 

where d L ±x is the Lebesgue measure on m + L 1 - induced from M d . 

Since dQ D A = 0, each subspace m + L 1 - for m G A either intersects the interior 
of P or is at least some distance e = e(P, L) > away from P. Hence we may 
replace 4> by the indicator [P] in (2.2.1). 

Recall from Section 1.2 that measuring volumes in m + L- 1 , we scale the volume 
form in L 1 - induced from W 1 so that the determinant of the lattice L 1 - fl Z d is 1. 
One can observe that det A provides the required normalization factor, so 

(det A) f [P] (x) d L ± (x) = vol d - k (P n (m + L L ) ) . 

The proof now follows. □ 

Remark. If dQ fl A 7^ the limit still exists but then for m G dQ n A the volume 
voLj-jt [P fl (m + -k" 1 )) is counted with the weight defined as follows: we find the 
minimal (under inclusion) face F of P such that m + L 1 - is contained in aff (F) and 
let the weight equal to the solid angle of P at F. 
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3. Exponential valuations 



Let V be a <i-dimensional Euclidean space, let A C V be a lattice and let A* be 
the reciprocal lattice. Let us choose a vector / G A* and let us consider the integral 



®i(P) = J exp{27rz(/, x)} dx, 



where dx is the Lebesgue measure in V. Note that for / = we have <&i(P) = 
<$> (P) = vol P. We have 

$ ; (P + a) = ex.p{2ni(l,a)}^i(P) for all a G V. 

It follows that <&i is a A-valuation on rational polytopes P G V. 

If I 7^ then the following lemma (essentially, Stokes' formula) shows that $j 
can be expressed as a linear combination of exponential valuations on the facets of 
P. The proof can be found, for example, in [B93b]. 

(3.1) Lemma. Let P G V be a full- dimensional polytope. For a facet Y of P , let 
drx be the Lebesgue measure on aff(r), and let pr be the unit outer normal to V. 
Then, for every I G V \ 0, we have 

J exp{27rz (/,£)} dx = 2ni\\l\\ 2 / ex P i 2 ^^^ x )i d r x > 

where the sum is taken over all facets V of P. 

Let F C P be an i-dimensional face of P. Recall that by lin(F) we denote the 
i-dimensional subspace of W 1 that is parallel to the affine hull aff(F) of F. We need 
the following result. 

(3.2) Theorem. Let P C V be a rational full- dimensional polytope and let t be a 
positive integer such that tP is a lattice polytope. Let e > be a rational number 
and let a G V be a vector. Let us choose I G A*. Then there exist functions 
fi(P, e, a, I; ■) : N — > C, i = 0, . . . , d, such that 

d 

(3.2.1) <f>i((n + e)P + a) = ^2f i (P,e,a,l;n)n i for all n G N 

i=0 

and 

(3.2.2) fi(P,e,a,l;n + t) = fi(P,e,a,l;n) for all neN 
and i = 0, . . . , d. 

Suppose that /d_fc(P, e, a, /; n) 7^ for some n. Then there exists a (d — k)- 
dimensional face F of P such that I is orthogonal to lin(F). 



Proof. Since 



$z(P + a) = exp {2ttz(/, a)} 



without loss of generality we assume that a = 0. We will denote fi(P, e, 0, /; n) just 
by /i(P,e,/;n). 

We proceed by induction on <i. For d = the statement of the theorem obviously 
holds. Suppose that d > 1. If / = then <E»z((to + e)P) = (n + e) d vo\P and the 
statement holds as well. 

Suppose that / ^ 0. For a facet V of P, let Ar = A n lin(r) and let Zr be the 
orthogonal projection of / onto lin(r). Thus A r is a lattice in the (d— 1) -dimensional 
Euclidean space lin(r) and l-p G Af, so we can define valuations on lin(r). Since 
tP is a lattice polytope, for every facet Y there is a vector «r G F such that 

lin(r) = aff(tr) - tw r and tw r e A. 

Let r' = T — ur, so r" C lin(r) is a Ar-rational (d — 1) -dimensional polytope 
such that tV is a Ar-polytope. We have 

(n + e)r = (n + e)r' + (to + e)u r . 

Applying Lemma 3.1 to (to + e)P, we get 

^ ((to + e)P) = £ V(r, /; n)$ Zr ((to + e)r') , 
r 

where 

V>(r, /; ra) = 2^pp ex P{ 2 ^(^ + e ) (*» w r> } 

and the sum is taken over all facets V of P. 
Since tur G A and / G A* , we have 

V>(I\Z;n + t) = tp(T,l;n) for all to G N. 

Hence, applying the induction hypothesis, we may write 

/i(P,e,/;n) = J]V(r,/;n)/ i (r / ,e,/r;n) for all togN 
r 

and i = 0, . . . , d - 1 and f d (P,e,l;n) = 0. Hence (3.2.1)-(3.2.2) follows by the 
induction hypothesis. 

If fd-k(P, e 7 /; n) 7^ then there is a facet T of P such that fd-k(X\ ^rS ^) ^ 0. 
By the induction hypothesis, there is a face F' of T' such that dimP' = d — k, and 
lr is orthogonal to lin(P'). Then P = P' + w r is a (<i — /c)-dimensional face of P, 
lin(P') = lin(P) and / is orthogonal to lin(P), which completes the proof. □ 
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4. Proof of Theorem 1.3 



First, we discuss some ideas relevant to the proof. 

(4.1) Shifting a valuation by a polytope. Let V be a <i-dimensional real vector 
space, let A C V be a lattice, and let v be a A- valuation on rational polytopes. Let 
us fix a rational polytope R C V. McMullen [Mc78] observed that the function \i 
defined by 

fJL(P) = u(P + R) 

is a A-valuation on rational polytopes P. Here "+" stands for the Minkowski sum: 

P + R={x + y: xeP,yeR,y 

This result follows since the transformation P \ — > P + R preserves linear depen- 
dencies among indicators of polyhedra, cf. [MS83] . 

Let t be a positive integer such that tP is a lattice polytope. McMullen [Mc78] 
deduced that there exist functions Vi(P, R; ■) : N — ► C, i = 0, . . . , d, such that 

d 

v (nP + R) = (P, R; n)n l for all tiGN 

and 

Vi (P, R; n + t) = Vi{P, R; n) for all tiGN. 

(4.2) Continuity properties of valuations E and E L . Let R C M, d be a full- 
dimensional rational polytope containing the origin in its interior. Then for every 
polytope Pcl d and every e > we have P C (P + eRj . We observe that 

\(P + eR) f]Z d \ = \Pf]Z d \, 

for all sufficiently small e > 0. If P is a rational polytope, the supporting affine 
hyperplanes of the facets of nP for n G N are split among finitely many translation 
classes modulo Z d . Therefore, there exists S = 5(P, R) > such that 

| (nP + eR)nZ d \ = \nP nZ d \ for all < e < 5 and all n G N. 

We also note that for every rational subspace L C M d , we have 

lim E L (P + eR) = E L {P). 

We will use the perturbation P i — > P + eR to push valuations E and El into a 
sufficiently generic position, so that we can apply Lemmas 2.1-2.2 without having 
to deal with various boundary effects. This is somewhat similar in spirit to the idea 
of [BS05]. 
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(4.3) Linear identities for quasi-polynomials. Let us fix positive integers t and 
d. Suppose that we have a possibly infinite family of quasi-polynomials pi : N — > C 
of the type 



Pi 



(n) = ^2pi(l;n)n l for all ti£N, 



i=0 



where functions pi(l; •) : N — > C, i = 0, . . . , d, satisfy 

Pi(l; n) = pi(l; n + t) for all ti£N. 

Suppose further that p : N — > C is yet another quasi-polynomial 

d 

p(n) = ^^pi(n)n l where Pi{n + t) = Pi(n) for all tiGN. 

i=0 

Finally, suppose that q(-) : M + — >■ C is a family of functions and that 
pin) = lim ci(r)pi(n) for all n£N 

T > + 00 ^— — ' 

I 

and that the series converges absolutely for every n e N and every r > 0. 
Then we claim that for i = 0, . . . , d we have 

pAn) = lim > c;(rk(I;n) for all n6N 

/ 

and that the series converges absolutely for every n G N and every r > 0. 

This follows since Pi(n), respectively pi(l;n), can be expressed as linear combi- 
nations of p(m), respectively pi(m), for m = n, n + t, ■ • • ,n+td with the coefficients 
depending on m, n, t, and d only. 

Now we are ready to prove Theorem 1.3. 

(4.4) Proof of Theorem 1.3. Let us fix a rational polytope P C M, d as defined in 
the statement of the theorem. For L e C let Pl C L be the orthogonal projection 
of .P onto L and let A/, C L be the orthogonal projection of Z d onto L. 
Let a G int P be a rational vector and let 

R = P -a. 

Hence R is a rational polytope containing the origin in its interior. Let Rl denote 
the orthogonal projection of R onto L. 

Since P is a rational polytope and £ is a finite set of rational subspaces, there 
exists 5 = 5(P, R) > such that for all < e < S and all n G N, we have 

(1) (nP + eR)nZ d = nPf]Z d and d (nP + eR) n Z d = for all n G N 
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and for all L G £, we have 

(2) (nP L + ei2 L )nA L = nP L nA L and d(nP L + eR L ) nA L = for all ti£N, 

cf. Section 4.2. Let us choose any rational < e < 5. 
Because of (1), we can write 

d 

(3) \{nP + eR)nZ d \= s ^ j e l {P]n)n i for all tiGN 
and by Lemma 2.1 we get 

(4) UnP + eR) f)Z d \ = lim V exp{-7r||Z|| 2 /r}$ z (nP + eR), 

1 1 T > + 0O ^ ' V ' 

iei d 

where are the exponential valuations of Section 3. Since is a Z d -valuation, 
by Section 4.1 there exist functions fi(P, e, I; •) : N — >■ C, z = 0, . . . , d, such that 

d 

(5) ®i(nP + eR) = e, /; n)n l for nGl 

i=0 

and 

(6) fi{P,e,l;n + t) = fi{P,e,l;n) for all tiGN. 
Moreover, we can write 

nP + eR = nP + e(P - a) = (n + e)P - ea. 

Therefore, by Theorem 3.2, for i < k we have fd-i{P, e, I; n) = unless / G L F for 
some face F of P with dimF = d — k. 

Therefore, combining (3)-(6) and Section 4.3, we obtain for all < i < k and 
all n G N 

e d -i(P;n)= lim V exp {-7r||/|| 2 /r} e, /; n) 

= lim V exp{-7r||/|| 2 /r}/ d - l ( J P,e,/;n), 

r >+oo ^ — * 

l€[j LeC (LnZd) 

since vectors / G Z d outside of subspaces L G £ contribute to the sum. Therefore, 
for < % < k and all n G N 

(7) e d . l (P;n)= lim V exp {-7r||/|| 2 /r} e, /; n) 

£e£ ieLnz d 
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On the other hand, because of (2), by Lemma 2.2 we get for all L G C and all n G N 

(8) E L (nP + eR)= lim V exp {-tt\\1\\ 2 /t} ^ t (nP + eR). 

Since E L are Z d - valuations, by Section 4.1 there exist functions 
ei(P, e, L; •) : N — > Q, i = 0, . . . , d, such that 

d 

(9) E L (nP + eR) =^Te i (P,e,L-,n)n i for all tiGN 

i=0 

and 

(10) ei(P,e,L;n + t) = ei(P,e,L;n) for all tiGN. 
Combining (5)-(6) and (8)-(10), by Section 4.3 we conclude 

ed-i{P, e, L; n) = lim N exp {— 7r||/|| 2 /rj f^-dP, e, /; n) for all n£N. 

ieLnz d 

Therefore, by (7), for < % < k we have 

(11) e d _i(P;n) = ^(L)e d _,(P,e,L;n) for all tiGN. 

Lec 

Since £x is a Z d -valuation, there exist functions ej(P, L; •) : N — ► Q, z = 0, . . . , d, 
such that 

d 

(12) E L (nP) = e t( p i L ; n K for a11 n eN 

i=0 

and 

ei(P, L; n + t) = ei(P, L; n) for all ti£l. 

Let us choose an m G N. Substituting n = m, m + t, . . . , m + dt in (12), we obtain 
ej(P, L; m) as a linear combination of E^nP) with coefficients depending on n, m, 
t, and <i only. Similarly, substituting n = m,m + t, . . . ,m+ dt in (9), we obtain 
ei(P, e, L; m) as the same linear combination of E^{nP + eRj. Since volumes are 
continuous functions, in view of (2) (see also Section 4.2), we get 

lim El (nP + eR) = El(uP) for n = m,m + t, . . . , m + dt. 

Therefore, 

lim €i(P, e, L; m) = ei(P, L; m) for all m G N. 

e >0+ 

Taking the limit as e — ► +0 in (11), we obtain for < i < k 

e d _ i (P;n) = J2^(L)e d _ i (P,L;n) for all tiGN. 
Lec 

To complete the proof, we note that 

Vd-i{P, L; n) = v(L)e d _i(P, L; n). 
Lec 



□ 
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5. Summing up a polynomial over 
integer points in a rational polytope 

Let us fix a positive integer k and let us consider the following situation. Let 
Q C M fc be a rational polytope, let int Q be the relative interior of Q and let 
/ : ~§L k — > R be a polynomial with rational coefficients. We want to compute the 
value 

(5.1) Yl /( m )- 

meintQnZ k 

We claim that as soon as the dimension k of the polytope Q is fixed, there is a 
polynomial time algorithm to do that. We assume that the polytope Q is a given 
by the list of its vertices and the polynomial / is given by the list its coefficients. 
For an integer point m = (/ti, . . . , /Ufc), let 

x m = x?---x% k for x= (xi,... ,x k ) 

be the Laurent monomial in k variables x = (x±, . . . , Xk). We use the following 
result [BP99]. 

(5.2) The short rational function algorithm. Let us fix k. There is a poly- 
nomial time algorithm, which, given a rational polytope Q C M fc computes the 
generating function (Laurent polynomial) 

S(Q;x)= Yl x?n 

m£ int QnZ k 

in the form 

S(Q; x) = g ei (1 _ xbil) ... (1 _ x6ifc) , 

where £ Z h , bij E Z h \ {0} and q G Q. In particular, the number |/| of fractions 
is bounded by a polynomial in the input size of Q. 

Our first step is computing the generating function 

5(Q,/;x)= fi™)*™- 

meQnz h 

Our approach is similar to that of [D+04]. 

(5.3) The algorithm for computing S(Q,f;x). We observe that 

S(Q,/;x) = / - ' x «^) ^ x )- 

14 



We compute S(Q; x) as in Section 5.2. 

Let a = («i, . . . , ajt) be an integer vector, let bj = ■ ■ ■ ,Pjk) be non-zero 
integer vectors for j = 1, . . . , k and let 71, . . . , 7^ be positive integers. Then 



d 

<L j 



dxi ) (1 - x & i)ti •••(!- yi b k)ik 

ai (l-X b i)^---(l-X b 7fc + ^ 7 ^ Jl (l -x b 7j+1 n (1-X 6 -)T-' 



Consecutively applying the above formula and collecting similar fractions, we com- 
pute 

/ _d_ _d_\ 

; \ Xl d Xl '"' ' Xk dxJ (l-x b i)---(l-x b O 

as an expression of the type 

( 5 ' 3 ' 1 ) Z^i-xbi^-i ...(1- X b *)W 

where pj G Q, 7ji, . . . , 7jfc are non-negative integers satisfying 7^1 + . . . + -fjk < 
k + deg / and cij are vectors of the type 

cij = a + + . . . + /U fc 6 fc , 

where /ij are non-negative integers and \i\ + . . . + Hk < deg /. The number of terms 
in (5.3.1) is bounded by (deg f)°^ k \ which shows that for a k fixed in advance, the 
algorithm runs in polynomial time. 

Consequently, S(Q,f;x) is computed in polynomial time. 

Formally speaking, to compute the sum (5.1), we have to substitute x^ = 1 into 
the formula for S(Q, f;x). This, however, cannot be done in a straightforward 
way since x = (1, . . . , 1) is a pole of every fraction in the expression for S(Q, /; x). 
Nevertheless, the substitution can be done via efficient computation of the relevant 
residue of S(Q, /; x) as described in [B94a] and [BW03]. 

(5.4) The algorithm for computing the sum. The output of Algorithm 5.3 
represents S(Q, /;x) in the general form 

where €j G Q, a; G Z fc , bij G Z fc \{0}, and 7^ G N such that 7a + - • -+7ifc < /c+deg / 
for all z G I. 
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Let us choose a vector / e Q fc , / = (Ai, . . . , A&) such that (/, bij) 7^ for all i,j 
(such a vector can be computed in polynomial time, cf. [B94a]). For a complex r, 
let 

x(r) = (e rAl ,... ,e rAfe ) . 
We want to compute the limit 

hm G(r) for G(t) = 5(Q,/;x(t)). 

T >0 

In other words, we want to compute the constant term of the Laurent expansion of 
G(t) around r = 0. 

Let us consider a typical fraction 



(1 - X fe l)Tl • • • (1 - X^)^ 

Substituting x(r), we get the expression 



(5.4.1) 



(l_ e r/3i) 7 i . ..(l_ e T/? fc )7k' 



where a = (a, /) and = I) for i = 1, . . . , k. The order of the pole at r = 
is D = 7i + • • • + 7fc < + deg /. To compute the constant term of the Laurent 
expansion of (5.4.1) at r = 0, we do the following. 
We compute the polynomial 



i=0 



that is the truncation at r D of the Taylor series expansion of e ar . For % = 1, . . . , k 
we compute the polynomial pi(r) with degpi — D such that 

Pi(j) + terms of higher order in r 



1 - e T ^ 



at t = 0. Consecutively multiplying polynomials mod r D+1 we compute a poly- 
nomial u{t) with degtt = D such that 



qi^pj 1 (t) ■ ■ ■ p1 k (t) = u(t) mod r 



D+l 



The coefficient of r D in u(r) is the desired constant term of the Laurent expansion. 
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6. Computing E l (A) 

Let us fix a positive integer k. Let A C R d be a rational simplex given by the list 
of its vertices and let L C R d be a rational subspace given its basis and such that 
dimL = k. In this section, we describe a polynomial time algorithm for computing 
the value of £x(A) as defined in Section 1.2. 

Let : R d — ► L be the orthogonal projection. We compute the vertices of 
the polytope Q = pr(A) and a basis of the lattice A = pr(Z d ). For basic lattice 
algorithms see [Sc93] and [G+93]. 

As is known, as x G A varies, the function 

<p(x) = vol d - k (P x ) where P x = (A n (x + L^)) 

is a piece- wise polynomial on Q. Our first step consists of computing a decompo- 
sition 

(6.1) Q = \Ja 

i 

such that Ci C Q are rational polytopes (chambers) with pairwise disjoint interiors 
and polynomials fa : L — > R such that fa(x) = fax) for x E Ci. 

We observe that every vertex of P x is the intersection of x + L 1 - and some k- 
dimensional face F of A. 

For every face G of A with dimG = k — 1 and such that aff(G) is not parallel 
to L- 1 -, let us compute 

A G = {xeL: i + L^afffG) ^0}. 

Then is an affine hyperplane in L. Then number of different hyperplanes Aq is 
(p{k) anc | h ence i\y e y cu i q i n ^ a t most d°( h ) polyhedral chambers C^, cf. Section 
6.1 of [Ma02]. As long as a; stays within the relative interior of a chamber Cj, the 
strong combinatorial type of P x does not change (the facets of P x move parallel to 
themselves) and hence the restriction fa of onto Cj is a polynomial, cf. Section 
5.1 of [Sc93]. Since in the (d — k) -dimensional space x + L 1 - the polytope P x is 
defined by d linear inequalities, fa can be computed in polynomial time, see [GK94] 
and [B93b]. 

The decomposition 6.1 gives rise to the formula 

[Q] = EKM' 

where Qj are open faces of the chambers Ci (the number of such faces is bounded 
by a polynomial in d), cf. Section 6.1 of [Ma02]. Hence we have 

e l (a) = j2 E tw- 

j meQj(~\A 

We compute inner sums as described in Section 5. 
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7. Computing e d _ fc (A;n) 



Let us fix a an integer k > 0. We describe our algorithm, which, given a positive 
integer d > k, a rational simplex AcR d (defined, for example, by the list of its 
vertices), and a positive integer n, computes the number ed-k(A; n). 

We use Theorem 1.3. 

(7.1) Computing the set C of subspaces. We compute subspaces L and num- 
bers n(L) described in Theorem 1.3. Namely, for each (d — /c)-dimensional face F of 
A, we compute a basis of the subspace L F = (linF)- 1 . Hence dimL F < k. Clearly, 
the number of distinct subspaces L F is d°( k K We let C be the set consisting of the 
subspaces L F and all other subspaces obtained as intersections of L F . We compute 
£ in k (or fewer) steps. Initially, we let 

£ := |l f : F is a (d— /c)-dimensional face of a|. 

Then, on every step, we consider the previously constructed subspaces L e £, 
consider the pairwise intersections LnL F as F ranges over the (d — /c)-dimensional 
faces of A and add the obtained subspace L n L F to the set £ if it was not already 
there. If no new subspaces are obtained, we stop. Clearly, in the end of this 
process, we will obtain all subspaces L that are intersections of different L Fi . Since 
dimL F = k, each subspace L e C is an intersection of some k subspaces L Fi . Hence 
the process stops after k steps and the total number \C\ of subspaces is d°( k \ 

Having computed the subspaces L e £, we compute the numbers n(L) as follows. 

For each pair of subspaces Li,L2 G C such that L\ C L2, we compute the 
number /u(Li, L 2 ) recursively: if L\ = L 2 we let /(/(Li, L 2 ) = 1. Otherwise, we let 

j u(L 1 ,L 2 ) = - v(L u L). 



Lec 

L 1 CLCL 2 



In the end, for each L G £, we let 



/*(L)= ^ /*(L,Li). 



Hence /(/(Lj, Lj) are the values of the Mobius function on the set C partially ordered 
by inclusion, so 



Lec 



follows from the Mobius inversion formula, cf. Section 3.7 of [St97]. 

Now, for each L e C and m = n, n + t, . . . , n + td we compute the values of 
E^mA) as in Section 6, compute 

i/(toA) = n(L)E L (mA) 
Lec 

and find ^_fc(A;n) = ed-fc(A,n) by interpolation. 
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8. Possible extensions and further questions 

(8.1) Computing more general expressions. Let P C W 1 be a rational poly- 
tope, let a > be a rational number, and let u G M. d be a rational vector. One can 
show (cf. Section 4.1) that 

d 

((n + a)P + u) nZ d | = ^e,(P, a, u;n)n l for all n G N, 

i=0 

where ej(P, a, w; •) : N — >■ Q, i = 0, . . . , d, satisfy 

ej(P, a, it; n + £) = ej(P, a, m; n) for all n G N, 

provided t G N is a number such that £P is an integer polytope. As long as k is 
fixed in advance, for given a, u, n, and a rational simplex A C M d , one can compute 
ed_fc(A, a, tt; n) in polynomial time. 

(8.2) Computing the generating function. Let P C M. d be a rational polytope. 
Then, for every < i < d, the series 

+oo 

J2ei(P;n)t n 
n=i 

converges to a rational function fi(P;t) for \t\ < 1. 

It is not clear whether /d_fc(A;t) can be efficiently computed as a "closed form 
expression" in any meaningful sense, although it seems that by adjusting the meth- 
ods of Sections 5-7, for any given t such that \t\ < 1 one can compute the value of 
/d_/-(A; t) in polynomial time (again, k is assumed to be fixed in advance). 

(8.3) Extensions to other classes of polytopes. If k is fixed in advance, the 
coefficient ed-k(P', n ) can be computed in polynomial time, if the rational polytope 
P C M d is given by the list of its d + c vertices or the list of its d + c inequalities, 
where c is a constant fixed in advance. 

(8.4) Possible applications to integer programming. If P C ~§L d is a rational 
polytope given by the list of its defining inequalities, the problem of testing whether 
P fl Z d = is a typical problem of integer programming, see [G+93] and [Sc86] . 
Moreover, by a general construction of "aggregation" (see Section 16.6 of [Sc86]) 
the problem can be reduced in polynomial time for that for P = A. It would be 
interesting to find out whether efficient computation of ed-fc(A;n) can have any 
practical applications to testing whether A fl Z d = 0. 
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